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Abstract 

In this paper, a multiscale approach is introduced 
to classify the JERS-1 mosaic image over the 
Central African rainforest. A series of texture 
maps are generated from the 100 m mosaic image 
at various scales. Using a quadtree model and 
relating classes at each scale by a Markovian 
relationship, the multiscale images are classified 
from course to finer scale. The results are verified 
at various scales and the evolution of classification 
is monitored by calculating the error at each stage. 

Introduction 


The large scale mosaic of JERS-1 data over 
tropical forests during the Global Rain Forest 
Mapping (GRFM) project has provided the 
opportunity of applying a series of data analysis 
and classification schemes for mapping and 
monitoring land cover and land use (Saatchi et al., 
1999). The high resolution imagery and the large 
scale mosaicking allows the visualization of the 
landscape at both fine and coarse resolution. In 
other words, more than one scale of the landscape 
heterogeneity can be observed or analyzed. 
Therefore, the data and the image mosaic enhance 
the traditional photogrametry and remote sensing 
techniques where one resolution and thus one 
scale is often observed. In this paper, the Central 
Africa Mosaic is classified for general land cover 
types such as forest, degraded forest, nonforest, 
floodplain forest, ficxxlplain low vegetation, 
grassland savanna, woodland savanna, and water. 
In order to exploit the landscape multiscale 
features for extracting thematic information from 
the JERS-1 mosaic, we have developed a new 
supervised classifier. The classifier is based on a 
multiresolution segmentation method that exploits 
the first order texture features of SAR data at 


various scales. This method was initially 
developed for a nonsupervised segmentation 
approach (Bouman and Shapiro, 1994). In this 
study, we have modified this concept for a 
supervised approach by integrating the technique 
with a maximum likelihood estimation (MLE) 
classifier. The methodology is applied on the 
Africa mosaic. The classified map is compared 
with the AVHRR based maps at 1km resolution 
and the Landsat based maps and field data at 100 
m. At each scale, the accuracy of maps, 
differences and complementarities are analyzed. 

JERS-1 Africa Mosaic 

In late 1995, the JERS-1 satellite entered into its 
Global Rain Forest Mapping (GRFM) phase to 
collect high resolution SAR data over the entire 
tropical rainforest. SAR data over central and 
west Africa, from eastern coast of Kenya and 
Madagascar, to Liberia and Guinea in the west, 
were acquired from January to March, 1996. The 
area covered lies between 9° N and 9° S latitude, 
extends approximately 6000 km along the 
equator, and amounts to about 2200 NASDA 
processed scenes or 9 million km 1 . The Congo 
River basin, about 3.5 million km 2 , was also 
covered during the October and November 1996 
(the high water season of the river). Madagascar 
data were acquired in January' 1997. The entire 
data acquired over the African continent have been 
processed and calibrated. Figure 1 shows tire 
geographical area and the mosaic of JERS-1 SAR 
data during the GRFM project over Central and 
West Africa and tropical rain forest region. JERS- 
l SAR is an L-band spacebome SAR system 
launched by the National Space Development 
Agency of Japan (NASDA) in February, 1992- 
The system operates at 1.275 GHz with horizontal 
polarization for both transmission and reception. 
The spatial resolution of the system is 18 m in 


1 


both azimuth and range. The swath width is 75 
km and the incidence angle of radar at the center 
of swath is 38.5°. The single-look images have 
4.2 m pixel spacing in azimuth and 12.5 m in 
range and the standard three look image has 12.5 
m pixel spacing in both azimuth and ranjge. 
JERS-I covers the global land surface for several 
applications such as land survey, agriculture, 
forestry, fisheries, environmental protection, 
disaster prevention and coastal monitoring. The 
satellite flies on sun-synchronous orbits 568 km 
above the Earth surface with a recurrent period of 
44 days. 




Figure 1. Geographical region of JERS-1 data 
take during GRFM project and the SAR mosaic 
image. 


Approach 

The objective of using multiscale texture based 
classification is to segment an image into regions 
of distinct statistical behavior. An implicit 
assumption in this procedure is that the statistics 
of each region are stationary and that each region 
represents a landscape feature that extends over a 
significant area (Bouman and Shapiro, 1994). 
Unlike an unsupervised classification which 
allows the segmentation based inherent statistics 
of the image, the supervised approach forces the 
segmentation to statistics developed by training 
data sets created over desired land cover types. 

In a nutshell, the approach includes the following 
steps: 


(1) A sequence of images is produced from JERS- 
1 mosaic where each successive image has half 
the resolution of its predecessor. These images 
are texture measures such as mean, coefficient of 
variation, energy, entropy, variance, contrast, 
kurtosis, and skewness are developed by using 
2x2, 4x4, 8x8 windows over 100 meter mosaic 
image to create 200 m, 400 m, and 800 m images 
respectively. 

(2) A training data set for desired land cover types 
are created by using a 8x8 window over the 100 
meter resolution mosaic image. The choice of 8x8 
window is dictated by the resolution of the 
coarsest image and the fact that larger windows 
may introduced mixed pixel statistics in training 
data set. 

(3) At each scale a supervised maximum 
likelihood technique is used as a base classifier to 
classify the image to land cover types. This 
derived land cover map can be verified by using 
course vegetation maps available over the 
continent 

(4) The sequence of random fields or class labels 
from coarse to fine scale follow a Maikov chain. 
This relation is imbedded in the MLE classifier 
used at the finer scale which allows segmentation 
of finer scale images with a priori knowledge of 
courser resolution classified map. 

(5) Class movement and accuracy are monitored at 
each scale. 

To capture mathematically the scale-to- 
scale relationship between classes or image 
statistics, we use a hierarchical approach known 
as quadtree model in which a class at scale n is 
related to its corresponding class at scale n+1 by 
using a Markov chain approach. In other words, 
assuming that the finest resolution is defined as 

scale n=0, the probability distribution erf - a class <o 

from a set of classes Q at this scale can be 
obtained by the following relationship: 
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where l is scale index, and P(.) is the probability 
density function of a class. Given the above 



relationship for class distribution along scales, the 
MLE distance at each scale can be modified such 
that it would include the probability density 
functions at each scale: 
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wjjere gj (x) is the distance measured for each 

pixel x, I and j are the class types and p and £ are 
the mean and covariance vector and matrices 
developed for each class from the input texture 
images. Note that at 100 meter resolution or scale 
n=0, the input image to the MLE classifier is only 
the backscatter image. In this case, the above 
relation will take a scalar form. 





Figure 2. Illustration of multi resolution sequences 
of three images with pixel-to-pixel quadtree 
branches. 

Discussion 

One of the challenges of employing the MLE 
approach with the quadtree pyramid approach is to 
define the transitional probabilities of classes at 
each scale on the tree structure (the first term in 
equation 4). Since the transition probabilities are 
defined on the tree structure, each pixel class is 


only dependent on a single pixel at the coarser 
scale. To model the probabilities in equation (4), 
wc have used an autoregresive type transition 
function and estimated its coefficients by using a 
preliminary fine-to-coarse scale classification 
(Bouman and Shapiro, 1994; Irving et al., 1997). 
These log conditional probabilities are computed 
and stored in image form at each scale. The 
number of points in the pyramid is approximately 
4N/3 where N is the number of pixels at the finest 
scales. If we choose M as a fixed number erf 
classes for all scales, the total storage required is 
approximately 4MN/3. 

The procedure shows that the classification at each 
scale takes into account the information 
transferred to it by the classification erf the 
previous scale. This method is suitable for large 
scale JERS-1 mosaic of Central Africa 
(approximately 35000 x 25000 pixel image). At 
the course resolution, the image size is smaller 
and therefore the classification and validation is 
faster. The hierarchical nature of the classifier, 
insures that at the finer scale the derived land 
cover map is improved. At each scale the 
classification can be compared and verified with 
maps and field information and the resulting map 
is used as an initial condition at the next finer 
resolution until individual pixels at 100 meter 
mosaic are classified. In this study, we compare 
the results with AVHRR 1 km land cover map of 
Central Africa and discuss the accuracy at each 
stage of classification. 
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